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Abstract 

We report on results obtained with the Event Shape Engineering technique applied to Pb-Pb colli¬ 
sions at v / Snn = 2.76 TeV. By selecting events in the same centrality interval, but with very different 
average flow, different initial state conditions can be studied. We find the effect of the event-shape 
selection on the elliptic flow coefficient vi to be almost independent of transverse momentum pj, as 
expected if this effect is due to fluctuations in the initial geometry of the system. Charged hadron, 
pion, kaon, and proton transverse momentum distributions are found to be harder in events with 
higher-than-average elliptic flow, indicating an interplay between radial and elliptic flow. 
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1 Introduction 

Results from Lattice Quantum Chromo-Dynamics |Q] |2| predict the existence of a plasma of deconfined 
quarks and gluons, known as the “Quark Gluon Plasma” (QGP). This state of matter can be produced in 
the laboratory by colliding heavy nuclei at relativistic energies |3j[4}0. The QGP was found to behave as 
a nearly perfect liquid and its properties can be described using relativistic hydrodynamics (for a recent 
review, see ( 6 l). The current experimental heavy-ion programs at Brookhaven’s Relativistic Heavy Ion 
Collider (RHIC) and at CERN’s Large Hadron Collider (LHC) are aimed at a precise characterization of 
the QGP, in particular of its transport properties. 

The system created in a heavy-ion collision expands and hence cools down, ultimately undergoing a 
phase transition to a hadron gas, which then decouples to the free-streaming particles detected in the 
experiments A precision study of the QGP properties requires a detailed understanding of this ex¬ 
pansion process. If the initial geometry of the interaction region is not azimuthally symmetric, a hydro- 
dynamic evolution of a nearly-ideal liquid (i.e. with a small value of the shear viscosity over entropy ratio 
tj/s) gives rise to an azimuthally anisotropic distribution in momentum space for the produced particles. 
This anisotropy can be characterized in terms of the Lourier coefficients v„ of the particle azimuthal dis¬ 
tribution Q. The shape of the azimuthal distribution, and hence the values of these Lourier coefficients, 
depend on the initial conditions and on the expansion dynamics. The geometry of the initial state fluctu¬ 
ates event-by-event and measurements of the resulting v„ fluctuations pose stringent constraints on initial 
state models. A quantitative understanding of the initial geometry of the produced system is therefore 
of primary importance | 6 ). A number of different experimental measurements and techniques have been 
proposed to disentangle the effects of the initial conditions from QGP transport, including measurements 
of correlations of different harmonics (8), event-by-event flow fluctuations lf9l ITOl ITTl fl2l and studies in 
ultra-central collisions lfl~3l fT4ll . Recent results from pp and p-Pb collisions at the LHC, moreover, sug¬ 
gest that hydrodynamic models may be also applicable to small systems lfT5llT6llT7l 182191. This further 
highlights the importance of studying Pb-Pb collisions with more differential probes, to investigate the 
interplay between the initial conditions and the evolution, in the system where the hydrodynamic models 
are expected to be most applicable. 

One of the new tools for the study of the dynamics of heavy-ion collisions is the “Event Shape Engi¬ 
neering” (ESE) lf20l . This technique is based on the observation that the event-by-event variation of the 
anisotropic flow coefficient (v„) at fixed centrality is very large MT2l . Hydrodynamic calculations show 
that the response of the system to the initial spatial anisotropy is essentially linear for the second and third 
harmonic, meaning that the final state V 2 (and V 3 ) are very well correlated with the second (and third) 
order eccentricities in the initial state for small values of r] /s 0 12Jj [22 1. These observations suggest 
a possibility to select events in heavy-ion collisions based on the initial (geometrical) shape, providing 
new opportunities to study the dynamics of the system evolution and the role of the initial conditions. 

The ESE technique is proposed to study ensemble-averaged observables (such as V 2 and inclusive particle 
spectra) in a class of events corresponding to the same collision centrality, but different v n values. In 
this paper events are selected based on the magnitude of the second order reduced flow vector <72 (see 
Sect. 13. II) . The technique was recently applied to study correlations between different flow harmonics 
in the ATLAS experiment 1231 . In this paper we present the results on elliptic flow and charged particle 
specta in Pb-Pb collisions at ^Ann = 2.76 TeV obtained with ESE technique. The events selected with 
the ESE technique are characterized by the measurement of V 2 , to quantify the effect of the selection on 
the global properties of the event. In order to search for a connection between elliptic and radial flow 
the effect of the ESE selection on the inclusive transverse momentum distribution of charged hadrons, 
pions, kaons and protons is then studied. The results are presented for primary charged particles, defined 
as all prompt particles produced in the collision including all decay products, except those from weak 
decays of light flavor hadrons and of muons. The differential measurement described in this work could 
provide important constraints to identify the correct model for initial conditions and for the determination 
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of transport properties. The development of flow in hydrodynamical models is driven by the pressure 
gradients and anisotropy in the initial state. A correlation between anisotropic and radial flow may stem 
from the specific fluctuation pattern in the initial state and/or can be produced in the final state depending 
on the bulk and shear viscosity of the system ||7j]. 

A few important caveats, which can affect the selectivity of the ESE technique, have to be kept in mind 
in this study. First, the discriminating power of the q 2 selection depends on the multiplicity and \>2 value 
in the psudorapidity, rj, region where it is computed and on the intrinsic resolution of the detector used 
for the measurement. Second, non-flow effects (such as resonance decays, jets, etc. 1221 1 could bias the 
q 2 measurement. In this work we discuss both aspects in detail, making use of different detectors with 
different intrinsic resolution and different 7 ] coverage. 

The paper is organized as follows. In Sect. [2] a brief review of the ALICE detector and of the data 
sample is presented. In Sect.[3]the analysis technique, with an emphasis on the event selection and on the 
particle identification strategy, is discussed. The results are presented in Sect. |4] Their implication for 
the hydrodynamic interpretation is discussed in Sect. [5] Finally, we come to our conclusions in Sect. [6] 

2 ALICE detector and data sample 

The ALICE detector at the CERN LHC was designed to study mainly high-energy Pb-Pb collisions. It is 
composed of a central barrel (|fl| < 0.8 for full-length tracks), containing the main tracking and particle 
identification detectors, complemented by forward detectors for specific purposes (trigger, multiplicity 
measurement, centrality determination, muon tracking). A detailed description of the apparatus can be 
found in lt24l . The main detectors used for the analysis presented in this paper are discussed below. 

The main tracking devices in the central barrel are the Inner Tracking System (ITS) and the Time Projec¬ 
tion Chamber (TPC). They are immersed in a 0.5 T solenoidal field. The ITS is the detector closest to the 
interaction point. It is a six-layer silicon tracker with a very low material budget (~ 7% of one radiation 
length X( } ). The ITS provides information on the primary interaction vertex and is used to track particles 
close to the interaction point, with the first layer positioned at a radial distance of 3.9 cm from the interac¬ 
tion point and the sixth one at 43 cm. It can measure the transverse impact parameter (DCA rv ) of tracks 
with a resolution of about 300 (40) jum, for transverse momentum pt = 0.1 (4) GeV/c, allowing the 
contamination from secondary particles to be significantly reduced. The TPC lf25l is a large-volume gas 
detector (external diameter 5 m) which measures up to 159 space points per track, providing excellent 
tracking performance and momentum resolution (o PT / p\ ~ 6% at pj = 10 GeV /c) l26l . It is also used 
in this work to identify particles through the measurement of the specific energy loss, clE/ dr. The dE/dx, 
computed as a truncated mean utilizing only 60% of the available samples, has a resolution of ~ 5% in 
peripheral and ~ 6.5% in central collisions | [26l . At a radius of 3.7 m from the beam axis, the Time 
of Flight (TOF) detector measures the arrival time of particles with a total resolution of about 85 ps in 
Pb-Pb collisions, allowing a n/K (K/p) 2 a separation up to pj = 3 (5) GeV/c. The ALICE reconstruc¬ 
tion software performs tracking based either on the information from the TPC alone (TPC-only tracks) 
or on the combined information from the ITS and TPC (global tracks). The former have the advantage 
of an essentially flat azimuthal acceptance, and are used for V 2 and q 2 measurements. The latter provide 
better quality tracks (G in /pj ~ 1.5% at p-\ = 10 GeV/c) | [26l l. rejecting most of the secondary tracks. 
However, the acceptance and reconstruction efficiency of global tracks are not flat in azimuth and as a 
function of transverse momentum, mostly due to missing or inefficient regions of the ITS. These tracks 
are used for the p-\ distribution measurements. TPC-only tracks can be constrained to the primary vertex 
(reconstructed also using the ITS information) to provide better momentum resolution. 

The data used for this analysis were collected in 2010, during the first Pb-Pb run at the LHC, at a center- 
of-mass energy per nucleon ,/snn = 2.76 TeV. The hadronic interaction rate was of the order of 100 Hz, 
low enough to avoid any space charge distortion effects in the TPC (27|. The trigger was provided by the 
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VO detector Il28l . a pair of forward scintillator hodoscopes placed on either side of the interaction region, 
covering the pseudorapidity regions 2.8 < 77 < 5.1 (VOA) and —3.7 < T] < — 1.7 (VOC). Events were 
requested to have a signal in both sides of the VO, selecting roughly 0-90% most central collisions 1291 . 
The VO measures a signal whose average amplitude is proportional to the multiplicity of charged parti¬ 
cles. The VO acceptance times detection efficiency is approximately 90% and flat as a function of the 
particle pj, with only a small reduction to about 85% for p-\ < 300 MeV/c. Events are further selected 
offline using the timing information from the V0 and from a set of two forward Zero Degree Calorimeters 
(ZDCs), in order to reject contamination from beam-induced backgrounds (see ll29l. 1301 [3T 1 for a detailed 
discussion). After all selections, the event sample used in the analysis consists of about 16 million events. 


3 Analysis technique 

3.1 Centrality and the event shape selection 


The events which pass the basic selection described in Sect. [2] are divided in centrality classes based 
on the signal amplitude (proportional to the charged particle multiplicity) measured in the V0 detector, 
as described in |[29ll . Events in each centrality class are further subdivided into groups with different 
average elliptic event shapes based on the magnitude of the second order reduced flow vector q 2 If 22 l 
given as 


_ \Qi\ 
qi Vm 


(i) 


where M is the multiplicity and \Q 


2 — 


Q 2 


2,x 


+ Q 2 


2,y 


is the magnitude of the second order flow vector. 


In this paper, the flow vector Q 2 is calculated using the TPC or V0 detectors. In the TPC, tracks in the 
range 0.2 < pj < 20 GeV /c and 1 77 1 < 0.4 (to avoid an overlap with the rj region used for the V 2 and pj 
distribution measurements) are used to measure 


M M 

Q 2,* = Y, cos2( Pi • Qi,y = !>2<P„ 

i— 1 i= 1 


( 2 ) 


where <p, is the azimuthal angle of the i-th particle and M is the number of tracks in an event. 

In the forward rapidity region the V0 is used. This detector segmented into four rings, each consisting of 
8 azimuthal sectors, the flow vector is hence calculated as 


32 32 32 

Qi,x = Y vv !'COs2<P( , Q 2 , y = Y Vri-sin2<p,- , M = Y w h (3) 

i= 1 i= 1 i=l 


where the sum runs over all 32 channels, (p, is the angle of the center of the sector containing channel i, 
Wj is the amplitude measured in channel i and M is in this case the sum of the amplitudes measured in 
each channel. 

The discriminating power of q 2 depends on the magnitude of elliptic flow as well as on the track multi¬ 
plicity used in the q 2 calculation and on the performance of the detector, including the angular resolution 
or the linearity of the response to the charged particle multiplicity. The good resolution of the TPC and 
the large multiplicity at midrapidity are used to maximize the selectivity on q 2 . However, the ALICE 
central barrel acceptance enables only limited separation in pseudorapidity between the region used to 
calculate q 2 and the region used to calculate the observables (|A 77 1 = 0.1). This separation is introduced 
in order to suppress unwanted non-flow correlations, which typicaly involve only a few particles and 
are in general short-range. In order to further assess the contribution of non-flow correlations, the flow 
vector is also calculated using the V0 detectors. This leads to a separation of more than one unit in 
pseudorapidity between the two regions. 
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Fig. 1: (Color online) Distributions of q^ vc (top) and q)' i)C (bottom) as a function of centrality (left) and projections 
for two centrality classes, 0 - 1 % and 30 - 31 % (right). In each of the left panels the solid curve shows the average 
qi as a function of centrality, while the dashed and the dotted curves indicate the top 10% and the bottom 10%, 
respectively. 

In absence of correlations, the average length of Q_ grows as \/~M ll22ll : qi is introduced to remove this 
trivial part of the multiplicity dependence. In case of non-zero correlations (due to either collective flow 
or non-flow correlations), q 2 depends on multiplicity and on the strength of the flow as lf32l(22l 

<^)~ 1 + (( M - l ))<( v | + ^)) > ( 4 ) 

where the parameter & accounts for non-flow correlations, and the angular brackets denote the average 
over all events. 

In the case when the multiplicity is measured via the signal amplitude in the V0 detector, the first term 
in Eq. [4] (unity) has to be substituted by (e?)/(e,) 2 , where <?, is the energy deposition of a single particle 
i. The fluctuations in e, lead to an increase in the flow vector length and reduce the corresponding event 
plane resolution. 
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Fig. 2: (Color online) Effect of the c/) nc (c/{ l>( ) event shape selection on the q^ pc (c/) {|C ) distributions for events 
in the 30-31% centrality class. 

The q 2 distribution measured with the TPC (qV >C ) and VOC (c/Y 0C ) is shown in Fig. Q] as a function of 
centrality, and in two narrow centrality classes, 0-1% and 30-31%. As can be seen, q 2 reaches values 
twice as large as the mean value, as expected in case of large initial state fluctuations l20l . The q)' i)C 
is larger than qj pc , as the former is measured in a larger pseudorapidity window (integrating a larger 
multiplicity) and is sensitive to the fluctuations in <?,. Note also that the selectivity (discrimination power) 
of the two selection cuts is in principle different, due to the different detector resolution, and, in the case 
of VOC, smaller V 2 value at forward rj, fluctuations in e\ and large contribution of secondary particles. 

In the present analysis, the effect of the event shape engineering on V 2 and pj distributions is studied. 
The average flow and particle spectra are measured in the pseudorapidity range 0.5 < 1 77 1 < 0.8 in order 
to avoid overlap with the region used to calculate qj pc . The VOC selection is used to estimate the 
contribution of non-flow correlations to the event-shape selection, since it provides a large I] gap. As 
a further crosscheck, the analysis was also repeated using the V0A detector. The results obtained with 
V0A and VOC show a qualitative agreement with a better selectivity when the VOC is used (mostly due 
to the larger multiplicity in the acceptance of this detector and to the r/ dependence of the elliptic flow). 
We therefore report the results for events selected using qj pc and qj° c in this paper. 

Due to the limited statistics, the analysis has to be performed in relatively wide centrality classes (~ 10%). 
The length of q 2 changes within such large centrality intervals (Fig. 03- and a cut at a fixed value of c /2 
would introduce a dependence on the multiplicity that would obscure the effect of the event-shape se¬ 
lection. The <72 selection is therefore evaluated in narrow (1%-wide) centrality classes. The results 
presented in the next sections are obtained in two event-shape classes, corresponding to the 10 % of the 
events having the top (bottom) value of the <72 (estimated in the narrow centrality classes). In the fol¬ 
lowing, we refer to these two classes as “large—< 72 ” (90-100%) and “small—< 72 ” (0-10%) or, generically, 
as ESE-selected events. Conversely, we refer to the totality of data within a given centrality class as the 
“unbiased” sample. 

The correlation between < 72 pc and q^ c is illustrated for events in the 30-31% centrality class in Fig. [2] 
The left (right) panel shows the distribution of <72 measured with the TPC (VOC) for all events and for 
events in the large— q 2 and small—c /2 classes, selected with the VOC (TPC). The average c /2 changes by 
about 18% and 14% in the large —^2 and small—c /2 samples respectively. In order to control the effect 
of fluctuations in a given detector the detailed comparison of the results obtained with qj pc and q\ {)C is 
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crucial, as discussed in detail below. In order to disentangle the effect of the rj gap and of the <72 cut, the 
selection on qV >( is also adjusted such that the average flow measured at mid-rapidity is similar to the 
one in the large —<72 sample (Sect. [4]). 

The ESE becomes less selective in peripheral events regardless of the detector used to compute q 2 , due 
to the low multiplicity. This limits the present analysis to the 60% most central events. 

Space charge distortion effects in the TPC, which accumulate over many events could, in principle, bias 
the q 2 selection. In order to check for this and other possible instrumental effects it was verified that the 
results are not sensitive to the instantaneous luminosity. 


3.2 Elliptic flow measurement 


The elliptic flow, V 2 , is measured in the pseudorapidity range 0.5 < | t]| < 0.8 using the Scalar Product 
(SP) method If22l . according to: 


-I*-:- 

J Q 2 Q 2 * / M a M b 


(5) 


where M 2 ,(t=exp(i 2 <pyi) is the particle’s unit flow vector, <p& is the azimuthal angle of the A-th particle of 
interest, Q 2 is the flow vector and M is the multiplicity. The full event is divided in two independent 
sub-events, labeled as A and B, covering two different pseudorapidity ranges, 0.5 < t] < 0.8 and —0.8 < 
77 < —0.5. The particle’s unit flow vector M 21 is evaluated in the sub-event A while the flow vector Qi and 
the multiplcity M in the sub-event B and vice-versa, ensuring a pseudorapidity gap of | Arj | > 1 between 
the particle of interest and the reference charged particles, which suppresses the non-flow contribution in 
the calculation of V 2 {SP}. A flat acceptance in azimuth is achieved in this analysis selecting TPC-only 
tracks, constrained to the primary vertex. Tracks are required to have at least 70 clusters and a (% 2 ) < 4 
per TPC cluster (two degrees of freedom). Tracks with a transverse distance of closest approach to the 
vertex (computed before constraining tracks to the primary vertex) DCA ty > 2.4 cm or a longitudinal 
distance of closest approach DC A- > 3.2 cm are rejected to reduce the contamination from secondary 
tracks. The effect of secondary particles is corrected applying the same analysis procedure to Monte 
Carlo events, simulated with the AMPT event generator Il33l and propagated through a GEANT3 Il34l 
model of the detector. The V 2 {SP} computed using reconstructed tracks is then compared with the one 
computed with generated primary particles, and the difference (< 5%) is used as a correction factor. 


The uncertainty on the tracking efficiency was assessed with different track samples and selections: 
using a set of hybrid tracks, built from a combination of global and TPC only tracks to obtain a uniform 
azimuthal acceptance |[35l . using TPC only tracks not constrained to the primary vertex, varying the 
minimum number of TPC clusters required in the analysis from 70 to 50 (Track reconstruction in Tab. [Q 
and [ 2 ]) and weighting each track by the inverse of the (pr-dcpcndcnt) efficiency (Tracking efficiency). 


The procedure used to estimate the centrality percentiles leads to a ~ 1 % uncertainty in the definition of 
the centrality classes l29l . In order to propagate this uncertainty to the results presented in this paper, the 
measurement is repeated displacing the centrality percentile by 1%. For instance, the analysis in the 30- 
40% centrality class is repeated for the selection 30.3-40.4% (Centrality resolution). Moreover, tracks 
reconstructed at midrapidity (instead of the V0 signal) are used as the centrality estimator (Centrality 
estimator). 


The correction for the effect of secondary particles mentioned above is strongly model dependent, there¬ 
fore the difference between the V 2 estimated using generated AMPT particles and reconstructed tracks 
was used to estimate the corresponding systematic uncertainty, ~ 3.5% (0.7%) at p-\ = 0.2 (1.5) GeV/c 
( 1 Secondary particles). 

Moreover, the following systematic checks were considered: the dependence on the magnetic field con¬ 
figuration was studied analyzing separately samples of events collected with different polarities of the 
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Effect 

v 2 

v 2 large —^2 

v 2 small —<72 

Track reconstruction 

3.1% (0-20%) 

3.1% (0-20%) 

3.1% (0-20%) 


2.7% (20-60%) 

2.7% (20-60%) 

2.7% (20-60%) 


(p T =0.2 GeV/c) 

(p T =0.2 GeV/c) 

(p T =0.2 GeV/c) 


0.08% (0-20%) 

0.08% (0-20%) 

0.08% (0-20%) 


0.02% (20-60%) 

0.02% (20-60%) 

0.02% (20-60%) 


(p T = 1.5 GeV/c) 

(p T =l-5 GeV/c) 

(p T =l-5 GeV/c) 

Tracking efficiency 

0.07% 

0.35% 

0.14 % 

Centrality resolution 

0.21% 

0.35% 

0.35% 

Centrality estimator 

0.57% 

0.49% 

0.57% 

Secondary particles 

3.56% 

3.56% 

3.56% 


(p T = 0.2 GeV/c) 

(p T = 0.2 GeV/c) 

(p T = 0.2 GeV/c) 


0.8% 

0.8% 

0.8% 


(p T = 1.5 GeV/c) 

(p T = 1.5 GeV/c) 

(p T = 1.5 GeV/c) 

Magnetic field 

NS 

NS 

NS 

Charge 

NS 

NS 

NS 

Vertex 

NS 

NS 

NS 


Table 1: Summary of systematic errors on x’ 2 {SP} measurement. NS = not statistically significant. 


Effect 

v 2 large—^/unbiased 

v 2 small—^/unbiased 

Track reconstruction 

0.14% 

0.14% 

Tracking efficiency 

0.35 % 

0.21% 

Centrality resolution 

0.14% 

0.21% 

Centrality estimator 

0.14% 

0.07% 

Secondary particles 

0.07% 

0.35% 

Magnetic Field 

NS 

NS 

Charge 

NS 

NS 

Vertex 

NS 

NS 


Table 2: Summary of systematic errors on the v' 2 {SP} ratios. NS = not statistically significant. 

magnetic field {Magnetic field), analyzing positive and negative particles separately ( charge) and analyz¬ 
ing samples of tracks produced at different vertex positions: —10 < z vtx < 0 cm and 0 < z vtx < 10 cm 
(Vertex). These effects are found to be not significant. 

The systematic uncertainties in the vo measurements and in the ratios of V 2 in ESE-selected over unbiased 
events are summarized in Tab. Q] and [2] Only the checks and variations which are found to be statistically 
significant are considered in the systematic uncertainties If36l . Whenever the pj dependence of the 
uncertainty is not negligible, values for characteristic pj are given in the tables. 

3.3 Transverse momentum distribution measurement 

The measurement of the p-\ distributions uses global tracks, which provide good resolution on DCA rv 
(Sect. O), and hence good separation of primary and secondary particles. The track selection requires at 
least 70 clusters in the TPC and at least 2 points in the ITS, out of which at least one must be in the first 
two layers, to improve the DCA XV resolution. A /;-| -dcpcndcnt cut on the DCA rv , corresponding to 7 times 
the experimental resolution on DCA xy , is applied to reduce the contamination from secondary particles. 
Tracks with a yfi per point larger than 36 in the ITS and larger than 4 in the TPC are rejected. Finally, 
to further reduce the contamination from fake tracks, a consistency cut between the track parameters of 
TPC and global tracks was applied. For each reconstructed TPC track, the ^-difference between the 
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Effect 

N ch 

77* 

F 

p and p 

Track reconstruction 

< 0.035% 

0.07% 

0.07% 

0.07% 

Tracking efficiency 

0 .21% 

0 .21% 

0 .21% 

0 .21% 

Centrality resolution 

0.07% (pj > 1.5 GeV/c) 

0.07% (p T > 1.5 GeV/c) 

0.14% 

0.14% 

Centrality estimator 

0.35% 

0.35% 

0.35% 

0.35% 

PID 

- 

0.07% (pi > 1.5 GeV/c) 

0.07% 

0.07% 

Secondary particles 

< 0.035% 

< 0.035% 

< 0.035% 

0.07% 

Normalization 

1 .1% 

1 .1% 

1 .1% 

1 .1% 

Magnetic field 

NS 

NS 

NS 

NS 

Charge 

< 0.035% 

< 0.035% 

< 0.035% 

< 0.035% 

Vertex 

0.07% 

0.07% 

0.07% 

0.07% 

Table 3: Summary of systematic errors for the ratio of pj distributions between large -(/2 and unbiased events. NS 
= not statistically significant. 

Effect 

N c h 

n ± 

/F 

p and p 

Track reconstruction 

< 0.035% 

0.07% 

0.07% 

0.07% 

Tracking efficiency 

0.28% 

0.28% 

0.28% 

0.28% 

Centrality resolution 

0.07% (p T > 1.5 GeV/c) 

0.07% (pi > 1.5 GeV/c) 

0.14% 

0.14% 

Centrality estimator 

0.35% 

0.35% 

0.35% 

0.35% 

PID 

- 

0.07% (pi > 1.5 GeV/c) 

0.07% 

0.07% 

Secondary particles 

< 0.035% 

< 0.035% 

< 0.035% 

0.07% 

Normalization 

0 .6% 

0 .6% 

0 .6% 

0 .6% 

Magnetic field 

NS 

NS 

NS 

NS 

Charge 

< 0.035% 

< 0.035% 

< 0.035% 

< 0.035% 

Vertex 

0.07% 

0.07% 

0.07% 

0.07% 


Table 4: Summary of systematic errors for the ratio of pi distributions between small -^2 and unbiased events. NS 
= not statistically significant. 


track parameters computed using only the TPC information constrained to the vertex and the associated 
global track is required to be less than 36 lf37ll . Charged tracks are studied in the pseudorapidity window 
0.5 < 1 77 1 < 0.8, to avoid an overlap with the tff K calculation. 

Particles are identified using the specific energy loss dE/dx in the TPC and their arrival time in the TOF. 
The technique is similar to the one presented in 11311 . A track is identified as either a pion, a kaon or 
a proton based on the difference, in the detector resolution units, from the expected energy loss and/or 
time of flight n<Jp lD (with i being the particle identity under study). Below p-\ = 0.5 GeV/c, only the 
TPC information is used («(7p ID = /za| pc ). For larger pj, the TPC and TOF information is combined 

using a geometrical mean: nOp lD = \J (/za-' r p C ) 2 + (no\- iW ) 2 . Tracks are required to be within 3(7 pid of 
the expected value to be identified as n ± , K + or p (p). In the region where the 3(7 pid identification 
bands of two species overlap, the identity corresponding to the smaller ziOpid is assigned. This technique 
gives a good track-by-track identification in the following pi ranges: 0.2 < p-y < 4 GeV/c for n ± , 
0.3 < pi < 3.2 GeV /c for K + , 0.5 < pj < 4 GeV/c for p (p). The misidentification of tracks is below 
4% for pions, 25% for kaons and 10% for protons in those ranges. Further discussion on the ALICE 
Particle Identification (PID) performance can be found in ll26l 1381 . The results for identified particles 
are provided in the pseudorapidity range 0.5 < 1 77 1 <0.8. However, in the case of the q^ oc selection 
the results were also studied at mid-rapidity |y| < 0.5. Results for positive and negative particles are 
consistent. In the following, “pions”, “kaons” and “protons”, as well as the symbols “ 7 r”, “K” and “p”, 
refer to the sum of particles and antiparticles. 
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The results for the spectra in ESE-selected events are presented in terms of ratios between the distribu¬ 
tions measured in the large— q 2 (small— qf) samples and the unbiased sample. The unbiased spectra have 
already been reported in Il37ll38l . Most of the corrections (and uncertainties) cancel out in these ratios, 
allowing for a precise determination of the effect due to the event-shape selection, as discussed in detail 
below. The uncertainties can mostly arise due to effects that depend on the local track density, which are 
found to be small l39ll . 

The systematic uncertainties are summarized in Tab. [3] and Tab. [4] As mentioned before, only the checks 
and variations which are found to be statistically significant are considered in the systematic uncertain¬ 
ties If36l . 

The systematic uncertainty related to the tracking is estimated varying the track selection cuts. Instead of 
the standard TPC cluster cut, at least 120 (out of 159) pad-rows hits in the TPC and a fraction of shared 
clusters in the TPC < 0.4, are required ( Track reconstruction in Tab. [3] and [4]). 

The possible effect of a track-density-dependent efficiency (which would influence in a different way 
events with the large— and small— q 2 selection) is investigated using simulations based on the AMPT 
event generator Il33l and a parametric event generator tuned to reproduce the ALICE spectra and V 2 mea¬ 
surements lf39ll . This effect leads to an uncorrelated systematic error of about 0.2% and a normalization 
error of 0.4% ( Tracking efficiency). 

The uncertainty on the centrality is estimated varying the definitions of centrality classes by 1% and 
using tracks as the centrality estimator. These checks lead to an uncorrelated uncertainty of about 0.1% 
and 0.35%, respectively and a normalization uncertainty below 1% in the ratios of spectra ( Centrality 
resolution and Centrality estimator). 

The systematic effect related to the particle identification is studied performing several variations to the 
PID approach described above. The nO\>\\y cut is varied between 2 and 4. Alternatively, if a track is 
consistent with more than one particle assignment within the nO\>\\y cut, double counting is allowed. As 
compared to the standard strategy where only the identity closest to the measured nCpiD is selected, this 
approach leads to a slightly larger contamination from misidentffied tracks, but also to a larger efficiency. 
Finally, an exclusive «Cpid strategy was used, which drastically reduces misidentffication: a particle is 
only accepted if it is compatible with only one mass hypothesis at 3<7pid- As a further cross-check, a 
Bayesian approach l26l was also considered. This method allows for better control of contamination at 
high p \. Overall, the uncertainty related to the particle identification strategy is less than 0.1% {PID). 

The effect of secondary particles depends on the pj distribution of weakly-decaying primary particles, 
and could be different for the large— and small —qi samples. This effect is estimated to be at most 
~ 0.1% for protons with the TPC ESE selection and negligible in all other cases {Secondary particles). 

Possible effects related to the magnetic field and to the charge state are addressed studying separately 
events collected with different magnet polarities {Magnetic field) and different charges {Charge), as in the 
case of the V 2 {SP} measurement. Particles produced at different longitudinal position cross a different 
portion of the detector, with different reconstruction efficiency. The samples of events produced with a 
negative (—10 < z vtx < 0 cm) and positive (0 < z vtx <10 cm) longitudinal vertex coordinate with respect 
to the nominal interaction point were studied separately {Vertex). 

4 Results 

4.1 Charged particle elliptic flow 

The event-shape selection is studied in Fig. [3} where the V 2 {SP} as a function of p-\ is reported for the 
unbiased and ESE-selected samples, with both the (| 77 1 < 0.4) and qf oc (—3.7 < T] < —1.7) selec¬ 
tions in different centrality classes. Figure [4] shows the ratio between the V 2 measured with the large —72 
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Fig. 3: (Color online) Measurement of V 2 {SP} as a function of p-y in different centrality classes for the unbiased, 
the large —<72 and the small —<72 samples. Only statistical uncertainties are plotted (systematic uncertainties are 
smaller than the markers). 
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Fig. 4: (Color online) Ratio of i' 2 {SP} in the large —qi and small— q 2 samples to unbiased sample. Only statistical 
uncertainties are plotted (systematic uncertainties are smaller than the markers). 
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Fig. 5: (Color online) Comparison between the effect of the event-shape selection obtained with the standard VOC 
and with the tuned TPC selections (see text for details), in the centrality class 30-40%. Top: V 2 {SP}, bottom: 
ratios to the unbiased sample. Only statistical uncertainties are plotted (systematic uncertainties are smaller than 
the markers). 

(small— qi) selection and the unbiased sample. Selecting the 10% highest (lowest) qj pc samples leads 
to a change of 30-50% in the V 2 {SP} measured, depending on centrality. The change is smaller (~ 10- 
25%) in the case of q^ QC -based selection, as compared to the qJ FC case. As already indirectly inferred 
from the difference between 2 nd and 4 th order flow cumulants V 2 { 2 } and V 2 { 4 } in llT2ll . the elliptic flow 
response of the system to geometry fluctuations is almost independent of p-\ . For all centralities, the 
change observed in Fig. |4] depends indeed weakly on pj, up to at least 4-5 GeV/c. This indicates that a 
cut on <72 selects a global property of the event, likely related to the initial shape in the overlap region. 
The only exception to the previous observation is the 0-5% centrality class, where for the qj pc selection 
an increasing trend with pj is observed. In this centrality class the mean value of V 2 is small, due to 
the almost isotropic shape in the initial state. Moreover, relative flow fluctuations are large in central 
collisions, with a pj dependence similar to the one shown in Fig. l4l liT2ll . The analysis of the pj spectra 
presented in Section l4~2l gives additional insight into the trend observed in Figure 0] 

For pj >4-5 GeV/c, the ratio ESE-selected/unbiased V 2 {SP} increases for the large —<72 selection. 
This trend is more pronounced for the <y{ l>( selection and for the most central and the most peripheral 
classes. A fit with a constant over the full pj range yields % 2 per degree of freedom values in the range 
2-6 (depending on centrality) for the qV‘ C selection and < 2 for the q^ oc selection. Fitting the ranges 
pj < 5 GeV/c and pj > 5 GeV/c with two different constants indicates an increase for the large —<72 
selection of order 5% and 10% for the q% oc and </j pc selections, respectively. This difference could be 
due to a small non-flow-induced bias. At high pj the V 2 is believed to be determined by the path-length 
dependence of parton energy loss fl2l . 

The difference between the qJ FC and </^ oc can be due to the different selectivity (see Sect. 13.11) . but also to 
a different contribution of non-flow correlations between the <72 and the V 2 measurements. Replacing the 
q FFC selection with the qj° c one changes both non-flow and selectivity at the same time. To disentangle 
these two contributions, the selectivity of the </ dT>c selection was artificially reduced. This is achieved 
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Fig. 6: (Color online) Centrality dependence of the average V 2 {SP} variation in the large—and small— <72 sam¬ 
ples. 

either relaxing the selection itself or rejecting a random fraction of tracks for the computation of cff K , 
while still selecting 10% of the events. It is found that selecting the class 65-100% for the large— q2 
sample (0-55% for the small— qi sample) with c/'( pc , or alternatively rejecting 70% of the TPC tracks, 
leads to an average variation of the V 2 {SP} in the range 0.2 < pj < 4 GeV/c comparable to the one 
obtained with the standard 10% q^ (K selection. The results are shown in Fig. [5] for the centrality class 
30-40%. Not only is it possible to find a cut which leads to the same average variation in V 2 {SP}, 
but the pj dependence is very similar in both cases. Rejecting randomly 70% of the tracks changes 
the selectivity of qV‘ c without affecting non-flow correlations between the cff K selection and V 2 {SP} 
measurement (as the 77 gap is not varied). Also in this case, it is found that the effect of the <72 selection 
does not depend on pj. A similar result, with the same value of the relaxed cut or fraction of rejected 
tracks, is found for the centrality interval 10-50%. Moreover, as it will be discussed in the next section, 
the same relaxed selections lead to the same effect on the pj distributions. 

These checks demonstrate that the selectivity of the cut is the main reason for the difference between the 
TPC and V0C selections. Due to the large p gap, the non-flow contribution is expected to be negligible 
in the case of the q^ oc selection. The agreement observed in Fig. [^indicates that, in the centrality classes 
10-50%, this is also the case for the qV‘ C selection in the range p \ < 5 GeV/c, a transverse momentum 
region dominated by hydrodynamic effects f3~8l . It is worth noticing that the ATLAS Collaboration 
measured a modification of the elliptic flow of ~ 35%, nearly independent of pj up to ~ 12 GeV/c in 
the 20-30% centrality class, while measuring V 2 and <72 with a pseudorapidity gap of 0.7 units |[23l . The 
increasing trend in the centrality class 0-5% is also observed in ll23ll PI 

To study the centrality and the dependence of V 2 {SP} in ESE-selected event classes, we quantified 

1 See auxiliary figures available on the ATLAS Collaboration web page 
https://atlas.web.cern.ch/Atlas/GR0UPS/PHYSICS/PAPERS/HI0N-2014-03/ 
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Fig. 7: (Color online) Average V 2 {SP} variation as a function of the absolute (left) values and self-normalized 
(right) values of the q^ c and q ^'° c for several centrality classes. 
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Fig. 8 : (Color online) Ratio of the p\ distribution of charged hadrons in the large— <72 or small— qi sample to the 
unbiased sample (qT >c and qj pc selections) in different centrality classes. Only statistical uncertainties are plotted 
(systematic uncertainties are smaller than the markers). 

the average change for each centrality class fitting the ratios in the range 0.2 < pj < 4 GeV/c with a 
constant]. The centrality dependence of the average change in the large —qi and small— q 2 selection is 
reported in Fig. [6j The trend obtained with the qj pc and qJ QC selections is very similar, except for the 
most central class 0-5%, where the average is influenced by the non-flat trend seen in Fig. [4] This once 
again reinforces the conclusion that the non-flow contamination is small also in the TPC selection case 
for the bulk of particles. The relative importance of non-flow changes with centrality. A large non-flow 
bias would therefore introduce a centrality dependence in the relative trend between the q ] -, PC and q^ )C 
selections, which is not observed. The dependence of the V 2 {SP} variation on ql P< and qJ QC is shown 
for the centrality classes 5-10%, 30-40% and 50-60% in Fig. [7] The left panel shows the absolute qi 

2 The result of the fit is numerically equivalent to the direct computation of the integrated V 2 in the range 0.2 < pj < 
20 GeV/c. 
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values on the x axis, while the right panel depicts the self-normalized values, defined as the average r /2 
value in ESE-selected events over the average qi values for all events in a given centrality class. The 
VOC selection spans a larger range but the TPC is more selective, as is clearly seen from the different 
slope of the TPC and VOC curves. In both cases the average q 2 reaches values twice as large compared 
to those in the unbiased sample, (Fig. [7] right). 

In summary, the observations reported in this section indicate that the ESE selects a global property of 
the collisions, as suggested by the flat modification in the V 2 as a function of pj. The q^ pc leads to a 
change twice as large than the corresponding q)' i)C selection. The difference between the two seems to 
be mostly due to the different discriminating power rather than to non-flow effects. 



Fig. 9: (Color online) Ratio of the pj distribution of charged hadrons in the large —qi or small —qi sample to the 
unbiased sample. Comparison between the effect of the event-shape selection obtained with the standard VOC and 
with the tuned TPC selections (see text for details), in the centrality class 30-40%. 


4.2 Transverse momentum distributions 

In order to study the interplay between the initial configuration of the system and the dynamics of the 
expansion of the fireball, the effect of the ESE selection on the single particle pj distribution is reported 
in Fig. [8] for the qV >( and q^ oc selections. As discussed in Sect. 13. li the reduced flow vector is calculated 
in the TPC detector in the pseudorapidity range 1 77 1 < 0.4. In order to avoid overlap between the qV’ ( and 
pj distribution measurements, only the region 0.5 < | t] | < 0.8 is used to measure the p\ distributions. 
This ensures at least 0.1 units of pseudorapidity separation between the qi and spectra measurements, 
thus suppressing the effect of short-range correlations. For consistency with the TPC analysis, the same 
pseudorapidity range is used in the case of the VOC selection. In the q^ oc case, it is also possible to study 
the spectra at mid-rapidity 1 77 1 <0.8 without any overlap with the qi measurement. The results agree 
within uncertainty with those in 0.5 < 1 77 1 < 0.8. 

The spectra in the large —<72 sample are harder than those in the small —<72 one. The ratio to the unbiased 
spectra reaches a maximum around pj = 4 GeV/c, and then stays approximately constant within large 
uncertainties. 

The effect of the selection is more pronounced in semi-central events (~ 30-50%), and decreases both 
towards more central and more peripheral collisions. This can be understood as due to the fact that the q 2 
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Fig. 10: (Color online) Ratio of the pj distribution of identified charged hadrons in the large —<72 sample to the 
unbiased sample for the f/| l>c (top) and c/Y 0C (bottom) selections. 


spans a larger dynamic range in semi-central collisions (Fig. [Hand Fig. [7}. In the most peripheral central¬ 
ity class studied in this paper (50-60%) the effect of the TPC-based selection is still very pronounced, 
while the q^ oc selection is less effective. This may indicate a small contamination from non-flow ef¬ 
fects in the most peripheral bin, consistent with observations discussed for the V 2 {SP} measurement in 
Sect. 14.11 In the most central class (0-5%) the modification of the spectrum is very small. This suggests 
that the trend observed in the same centrality class in Fig.[4]is likely to be dominated by flow fluctuations 
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Fig. 11: (Color online) Ratio of the pj distribution of identified charged hadrons in the small— q 2 sample to the 
unbiased sample for the qj pc (top) and r/) 0C (bottom) selection. 


rather than non-flow contributions. 

As in the previous section, we disentangle the effect of non-flow and q 2 selectivity either relaxing the 
q pl>( selection or randomly rejecting a fraction of the tracks. The relaxed cut and the fraction of rejected 
tracks tuned to reproduce the V 2 variation in 0.2 < pj <4 GeV/c in Sect. 14. II are used. Figure [9] shows 
that these selections yield results compatible with the standard q^ oc selection. A similar result (with the 
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Fig. 12: (Color online) Ratio of the pj distribution of identified charged hadrons in the large —qi (top) and 
small— q2 (bottom) sample to the unbiased sample (qJ PC selection), in 30-40% centrality class. Lines: ratio of 
the blast-wave parametrizations (see text for details). 

same relaxed cuts or fraction of rejected tracks) is found for all centralities up to ~ 50%, after which 
non-flow effects seem to become relevant. 

As discussed in Sect. 14. 11 we conclude that the effect of non-flow is small and that the main factor driving 
these observations is the average vi at mid-rapidity. 

The modification on the spectra of identified n, K, and p is reported in Fig. |T0] and Fig. QT] for differ¬ 
ent centrality classes. The same pattern measured in the case of non-identified hadrons is observed. 
Moreover, a clear mass ordering is seen: the modification is more pronounced for heavier particles. Con¬ 
versely, the spectra in the small —qo sample are softer. In the case of the VOC selection the analysis was 
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also repeated in the region |y| < 0.5, yielding consistent results. 

These observations suggest that the spectra in the large— q 2 (small—< 72 ) sample are affected by a larger 
(smaller) radial flow push. This hypothesis was tested with a blast-wave ll40ll study. A ratio of two 
blast-wave functions was used to fit the spectra ratios shown in Fig.[l0]and Fig. [TT] The parameters were 
initially fixed to the values from j38j, where they were tuned to describe the inclusive spectra of pions, 
kaons and protons. Then, the (/3 t) parameter of the numerator function was allowed to change (while 
keeping the overall integral of the function constant). The fit was performed as in If38l in the transverse 
momentum ranges 0.5-1 GeV/c, 0.2-1.5 GeV/c, 0.3-3 GeV/c for n, K, p, respectively. The agreement 
with the data is good, also outside the range used to determine the parameters, up to pj ~ 3GeV/ c. The 
fits yield the following result for the difference A (/3t) between the (Jij) parameter of the numerator and 
denominator function: A(/3t) = (0.41 ±0.03)% (large— q^) and A(/3t) = (—0.22 ±0.03)% (small— qi) 
for the centrality class 30-40%, as shown in Fig. [12] 

5 Discussion 

In this paper the first application of the Event Shape Engineering (ESE) |[20l to the analysis of ALICE 
data was presented. 

The results on the vijSP} measurement suggest that the ESE technique selects a global property of the 
collision, likely related to the eccentricity in the initial state. The measurement of pj spectra indicates 
that events with larger eccentricity show an increased radial flow. A correlation between elliptic and 
radial flow could be introduced either at the initial stage, due to the specific fluctuation patterns in the 
energy deposition, or during the hydrodynamic evolution of the system, due to an interplay of bulk and 
shear viscosity f7|. 

A Glauber Monte Carlo simulation was performed to estimate the possible correlation between the ini¬ 
tial eccentricity and azimuthally-averaged pressure gradients. In the model, the multiplicity of charged 
particles in the acceptance of the V0 detector, used to determine the centrality classes, is computed 
following |[29l . A “number of ancestors” N ancestors is derived from the number of participant nucleons 
(Npan) and binary collisions {N co u) as 

Nancestors — Imparl ± (1 f)N C oU■ ( 6 ) 

Each ancestor is assumed to produce particles following a negative binomial distribution with parameters 
taken from lf29l . 

The participant density, defined following 1411142 .91 1431 as N part /S, is used as a proxy for the magnitude 
of the pressure gradients. The average cross-sectional area S and participant eccentricity e are computed 
as 


5 = 4 KGyf Gy = 471 a 2 a 2 - a 2 , (7) 

Gy ~ Gy j K 2 - °) 2 ) 2 +4 G^, 

£ = -4- \ = 1 ——j-5-(8) 

Gy ± Gy CT 2 ± Gy 

where 

°X = C* 2 ) - (x) 2 , <*y = (/> - (y) 2 , Oxy = {xy)-{x)(y). (9) 

The unprimed coordinates are given in the fixed laboratory coordinate frame. Primed coordinates, x' and 
/, are calculated in the so-called participant coordinate system, rotated with respect to the laboratory 
coordinate frame such that the minor symmetry axis of the participant nucleon distribution coincides 
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Fig. 13: (Color online) Participant density as a function of the participant eccentricity estimated in a Glauber 
Monte Carlo model for central (top) and semi-central (bottom) collisions. 


with the x' direction. The normalization of the area is chosen such that for a Gaussian distribution the 
average density coincides with N part / S. 

Two narrow centrality classes, selected based on the simulated charged particle multiplicity, roughly cor¬ 
responding to 0-2% (central) and 30-32% (semi-central), are studied in Fig. [[3] The observed correlation 
between the density and the participant eccentricity is reminiscent of the correlation between radial flow 
and event shape measured in this paper. The average density in events with the 10% largest £ is about 1% 
(7%) larger than in events with the smallest £ for central (semi-central) collisions, qualitatively consistent 
with what is observed in Fig. [TO] and Fig. |TT] where the effect of the ESE selection is much stronger for 
semi-central collisions. This reinforces our conclusion that ESE is an effective tool to select the initial 
shape and density, thereby opening the possibility of further studies. 

A quantitative comparison would require a full hydrodynamical calculation. The correlation can in fact 
be modified by the transport in the hydrodynamic phase. In particular, it was shown [j7] l44l that in a 
system with a finite shear viscosity the flow coefficients, obtained for a given set of initial eccentricities, 
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Fig. 14: (Color online) Measurement of i' 2 {SP} as a function of pj, for the unbiased sample and for the large— q 2 
sample (top panel) and ratio between the large —<72 result over the unbiased result (bottom panel). Data points 
(full markers) are compared with AMPT Monte Carlo model (bands). Only statistical uncertainties are plotted 
(systematic uncertainties are smaller than the markers). 

are reduced as compared to the ideal hydrodynamics case. At the same time, shear viscosity increases 
the radial flow. In principle, bulk viscosity reduces the radial flow, reducing the correlation observed 
in this paper, but the latter effect was estimated to be negligible 04!. Therefore, the measurement we 
present in this paper is sensitive to the interplay of initial conditions and transport coefficients in the 
hydrodynamic phase. As such, it poses stringent constraints on hydrodynamic calculations, and it could 
allow the extraction of the value of average shear viscosity at the LHC. 

A study of the relation of the fluctuation in the initial size to the spectra was performed in Il45l 1461 
with a full hydrodynamic simulation. It was shown that the event-by-event fluctuations in the Glauber 
initial conditions lead to fluctuations in the initial size of the system that reflect in fluctuations of the 
radial flow and hence (pr). It is found that the relative (pj) fluctuations computed with Glauber initial 
conditions overestimate the data, indicating a strong sensitivity of event-by-event measurements on the 
initial conditions model. It is also shown that the (pj) fluctuations are not sensitive to the shear viscosity. 
The study in l45l l46l (fluctuations in (pr)), however, does not address the relation between the elliptic 
and radial flow. It may be expected that the present measurement will also be sensitive to the transport 
coefficient of the medium. 

In a recent series of theoretical studies B71 (48l l49ll . it was suggested to use the Principal Component 
Analysis to study flow fluctuations. It was argued that most of the current methods to study flow do 
not fully capture the complexity of the initial state. Indeed, the PCA studies revealed the presence of 
sub-leading flow components (arising from radial geometry excitations), which break the factorization 
of flow harmonics 0711481 . In particular, in |[49l it is argued that the sub-leading component of V 2 reflects 
a non linear mixing with radial flow, which could address the same physics as reported in this paper. 

A study of the relation of the fluctuation in the initial size to the spectra was performed in ll45l [46l 
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with a full hydrodynamic simulation. It was shown that the event-by-event fluctuations in the Glauber 
initial conditions lead to fluctuations in the initial size of the system that reflect in fluctuations of the 
radial flow and hence (pr). It is found that the relative (pj) fluctuations computed with Glauber initial 
conditions overestimate the data, indicating a strong sensitivity of event-by-event measurements on the 
initial conditions model. It is also shown that the (pj) fluctuations are not sensitive to the shear viscosity. 
The study in Il45ll46l (fluctuations in (pr)), however, does not address the relation between the elliptic 
and radial flow. It may be expected that the present measurement will also be sensitive to the transport 
coefficient of the medium. 

In a recent series of theoretical studies H71 1481 1491 . it was suggested to use the Principal Component 
Analysis to study flow fluctuations. It was argued that most of the current methods to study flow do 
not fully capture the complexity of the initial state. Indeed, the PCA studies revealed the presence of 
sub-leading flow components (arising from radial geometry excitations), which break the factorization 
of flow harmonics |[47ll48ll . In particular, in 1491 it is argued that the sub-leading component of V 2 reflects 
a non linear mixing with radial flow, which could be related to the same underlying physics phenomena 
reported in this paper. 

To further understand the observed effect, we studied it in AMPT, a model known to reproduce many 
of the flow observables measured at the LHC lf33ll . This model is based on HIJING lf50l to describe 
the initial conditions and on the Zhang’s parton cascade 1511 to describe the partonic evolution. The 
string melting configuration, described in If52il . is used. To assess the impact of the detector resolution 
on the q 2 selection, the simulated AMPT events were transported through the ALICE apparatus using 
the GEANT f34ll transport model. The qi was computed using either the reconstructed Monte Carlo 
tracks (q'^) or the generated primary particles in the same kinematic range (<yf cn ). The elliptic flow and 
the transverse momentum distribution are calculated using generated Monte Carlo particles. Since the 
charged particle multiplicity distribution is different in AMPT and data, the q 2 selection is calibrated 
in the model as a function of multiplicity. The results are shown in Fig. [IT] for the charged hadrons 
elliptic flow and in Fig. [15] for the transverse momentum distribution of charged hadrons. Using either 
qf c or q^ en does not introduce any significant difference on the effect of the selection. This indicates 
that detector resolution effects are negligible for the qj pc selection. The V0 detectors, on the other hand, 
have a coarser azimuthal resolution and are sensitive to fluctuations in the energy deposition of incident 
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particles. However, the study with the relaxed TPC selection discussed in Sect. [4] demonstrates that the 
properties of the ESE selected events are mostly determined by the average V 2 {SP} value. It is therefore 
advised that in any comparison of this data to theoretical models the selection in the model is tuned as to 
reproduce the average change in V 2 {SP} at mid-rapidity. 

The pj dependence of the elliptic flow observed in data is not reproduced in AMPT (top panel). This 
model reproduces however the magnitude of the modification as well as the flatness of the ratio as a 
function of pj. 

The effect of the ESE selection on the pj distribution of charged particles is well reproduced by AMPT 
below pj = 2 GeV/c, as shown in Fig. [15] However, the magnitude of the effect at intermediate pj 
(2 < pj < 6 GeV /c) is underestimated in AMPT. As previously observed for the V 2 measurement, a good 
agreement is observed between the selection based on g| en and q'^ . 

6 Conclusions 

In summary, the first application of the Event Shape Engineering (ESE) technique to Pb-Pb collisions 
data measured by ALICE at ^/snn = 2.76 TeV has been presented. 

The elliptic flow at mid-rapidity is observed to increase as a function of the <72 calculated in the central 
or forward rapidity regions. The modification of the V 2 coefficient as a function of pj is nearly flat 
below pj = 4 GeV/c, suggesting that this technique allows the selection of a global property of the 
collision, likely related with the geometry of the participant nucleons in the initial state. In the region 
above pj > 5 GeV/c a small increase is observed within the large statistical uncertainties, possibly due 
to a small non-flow contamination. In this transverse momentum range the elliptic flow is believed to be 
driven by the different path-length traversed in- and out-of-plane by high-px partons in the deconfined 
medium, rather than by the hydrodynamic evolution of the system. 

The pj -distributions of unidentified hadrons in the pj region (0 < pj < 5 GeV/c) are harder (softer) in 
event with large —qi (small— q^) values. 

Identified pions, kaons and protons show a similar behavior with a clear mass ordering in the ratio be¬ 
tween the large —<72 and the unbiased spectra, thus suggesting this effect to be due to a stronger radial 
flow in such events. Glauber Monte Carlo calculations reveal a correlation between the transverse par¬ 
ticipant density and the participant eccentricity which could be the origin of this effect. This indicates 
that at least part of the correlation is generated in the initial state. However, these measurements are also 
sensitive to the transport coefficients in the hydrodynamic evolution. A quantitative comparison would 
require a full hydrodynamic calculation and may provide stringent constraints both on shear and bulk 
viscosity. 
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